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Class Meeting #15: Classification of second order equations 



1. Review of Three Important Examples of PDEs 
Let's review some basic facts concerning the three PDEs we've examined in detail thus far. 



Equation 


Type 


Well-posed problems 


Features 


Au{x) = f{x) 


Elliptic 


Boundary value prob- 
lems: All of M" (with 
boundary conditions 
at oo); finite bound- 
aries under Dirichlet, 
Neumann, Robin, 
or Mixed boundary 
conditions 


mean value properties; 
maximum principle; Har- 
nack inequality 


dtu{t, x) — Au(t, x) = fit, x) 


Diffusive (parabolic): 


Initial value (Cauchy) 
problems: all of M" 
at t = 0; Initial + 
boundary value prob- 
lems: data at t = + 
Dirichlet, Neumann, 
Robin, or Mixed 
boundary conditions 


Infinite speeds of propaga- 
tion; smoothing properties; 
maximum principle, 
decay as t — )• oo for the 
global Cauchy problem 


—dfu{t, x) + Au{t, x) = f{t, x) 


Hyperbolic 


Initial value (Cauchy) 
problems: all of M" 
at t = 0; Initial + 
boundary value prob- 
lems: data at t = + 
Dirichlet, Neumann, 
Robin, or Mixed 
boundary conditions 


Finite speed of propagation; 
domain of dependence and 
influence; energy identities; 
order ^(i"")/^ decay as t ^ 
oo for the global Cauchy 
problem 



2. Motivating example 



Let's consider the following second-order linear PDE on M 



(2.0.1) 



Cu = A'^^dadpu + B"daU + Cu = 0. 



In ( |2.0.1 ), A,B,C are allowed to be functions of the coordinates 



X 



standard notation x° 
that A is symmetric: 



x"). We will also use the 



t. By the symmetry of the mixed partial derivatives, we can also assume 



2 
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(2.0.2) A^"" = A""". 

The question we would like to address at the moment is the following: what are the basic properties 
of solutions to (2.0.1)? Is this equation most like a Laplace, heat, or wave equation? That is, is 
(2.0.1) elliptic, diffusive, or hyperbolic? As we will see, the most important part of equation (2.0.1) 
in this context is the principal part A^^dad^u, which involves the top-order derivatives. 

To begin answering this question, let's start with a simple example on M^. Let's try to classify 
the following equation: 



(2.0.3) Cu = d^u - Adtd^u + 2dlu = 0. 

Note that it would be easy to answer our question if we were able to make a linear change of 
variables that eliminates the cross term —Adtd^u; the PDE would then look just like one of the 
other ones we have already studied. More precisely, let's try to eliminate the cross terms by making 
good choices for the constants a, b, c, d in the following linear change variables: 



(2.0.4a) 
(2.0.4b) 



t = at + bx, 
X = at + dx. 



In order to have a viable change of variables, we also need to achieve the following non-degeneracy 
condition from linear algebra: 



(2.0.5) 



ad — be ^ 0. 



(2.0.5) states the determinant of the above linear transformation is non-zero, and that the trans- 



formation is non-degenerate. 

Then using the chain rule, we have that 



(2.0.6a) 
(2.0.6b) 



dt 

dr. 



dt 

dt 
di 



dx 

r. dr+ —dx = adr+ cdx, 
dt ^ dt * 



dx 



Inserting (2.0.6a) - (2.0.6b) into (2.0.3), we compute that 



(2.0.7) 



Cu = (a^ - Aab + 2b')d~u + (2ac + Abd - Aad - Abc)d^dzu + (c^ - Acd + 2d^)d~u 



To make the cross term in (2.0.7) vanish, we now choose 



(2.0.8) 



0, 



d=l. 



Note that (2.0.8) also verifies the non-degeneracy condition (2.0.5). We remark that other choices 
would also have worked. In the new coordinates, we have that 



(2.0.9) 



Cu = d~u — 2d~u. 
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Dividing by —2, we see that the PDE (2.0.3) was actually a "standard" linear wave equation in 
disguise: 



ldlu + dl = 0. 



(2.0.10) 

Relative to the coordinates {t,x) , the "speed" associated to the wave equation (2.0.10) is a/2. 



Let's do another example. Consider the PDE 



(2.0.11) 



Cu = —2dfu — 2dtdxU — dlu + dxU = 0. 



Using (2.0.6a) - (2.0.6b) again, we compute that 



(2.0.12) Cu = {-2a' - 2ab - ¥)dfu + {-2ac - bd - 2ad - 2hc)d^xU + {-2c' - Acd - d')diu 



Choosing 

(2.0.13) 

we see that 



+ hdva + ddxU. 



1 

71' 



6 = 0, 



c = -1, 



d=h 



(2.0.14) 



Cu = —d-jrU — d^u + dxU. 



Thus, multiplying by —1, we see that (2.0.11) is really just a Laplace-like equation in disguise: 



(2.0.15) 



d~u + d~u — dxU = 0. 



Equation (2.0.11) is therefore elliptic. We remark that the first-order term in (2.0.15) does not 



affect the elliptic nature of the system. 

Let's do one final example. Consider the PDE 



(2.0.16) 



Cu =^ dfu — 2dtdxU + d'u + d^u = 0. 



Using (2.0.6a) - (2.0.6b) again, we compute that 



(2.0.17) Cu 



Choosing 



2ab + b')d~u + {2ac + 2bd - 2ad - 2bc)dtdsu + {c' - 2cd + d')d~u 



+ bdfU + ddxU. 



(2.0.18) 

we see that 



a = 1, 



6 = 0, 



(2.0.19) Cu = d^u-dxU. 

Thus, (2.0.16) is equivalent to 
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(2.0.20) -d^u + dfu = 0. 



Now observe that (2.0.20) is just the standard heat equation, with the variable x playing the role 



of "time" and t playing the role of "space." Equation (2.0.20) is therefore diffusive (parabolic). 



3. A GENERAL FRAMEWORK 

In this section, we will establish a general framework for classifying second order constant co- 
efficient scalar PDEs. The framework will cover the three examples from the previous section as 
special cases. The proof will reveal that the classification is intimately connected to the theory of 
quadratic forms from linear algebra. Throughout this section, we will use the notation 

(3.0.21) X = {x^,x\--- ,x''). 

As above, we will investigate PDEs of the form 

(3.0.22) Cu = A'^^d^dpu + B^'d^u + Cu = 0, 

where A'"'' = A"''. 

We begin by providing a simple version of Hadamard's classic definitions. 



Definition 3.0.1 (Hadamard's classification of second order scalar PDEs). Equation (3.0.22) 
is respectively said to be elliptic, hyperbolic, or parabolic according to the following conditions on 
the (1 + n) X (1 + n) symmetric matrix A : 

• All of the eigenvalues of A have the same sign - elliptic 

• n of the eigenvalues of A have the same (non-zero) sign, and the remaining one has the 
opposite (non-zero) sign - hyperbolic 

• n of the eigenvalues of A have the same (non-zero) sign, and the remaining one is - 
parabolic 

Remark 3.0.1. Many of the ideas in this section, including the definition above, can be generalized 
to include the case where A depends on (x), or even on the solution u itself; PDEs of the latter 
type are said to be quasilinear. 

We now state and prove the main classification theorem. 

Theorem 3.1 (Classification of second order constant-coefficient PDEs). Consider the 

following second order constant coefficient PDF 

(3.0.23) Cu{x) = A'"^dadi3u{x) + B"d^u{x) + Cu{x) = 0, 

where da ==^gf^- Then there exists a linear change of variables y'^ = MJ^x"" such that 



If all of the eigenvalues of A^'^ have the same (non-zero) sign, then (3.0.23) can be written 
as ±Cu = Ayu{y) + B^'^uiy) + Cu{y) = 0, where Ay = Y!l=o 

If n of the eigenvalues of A have the same (non-zero) sign, and the remaining one has the 



opposite (non-zero) sign, then (3.0.23) can be written as ±Cu = □j^n(y) + B°'-J^u{y) + 



Cu(y) = 0, where niy =^ (m ^)°'^ ^j^^ is the standard linear wave operator, and (m) ^ 
diag{—l, 1, 1, ■ ■ ■ ,1) is the standard Minkowskian matrix. 
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Ifn eigenvalues \^^\ ■ ■ ■ , A^"-* of A have the same (non-zero) sign, and the remaining one is 

= 0, then (|3X23|) can be written as ±Cu = B'^uiy^, y\-- - , y'')+Eti Sf^y""^ V^'-' ^ + 



X]r=i ~^'^{y^ -I -I ' ' ' + = 0. Furthermore, let v^^\v'^^\ ■ ■ ■ ,f he a corresponding 
diagonalizing unit-length co-vector basis. More precisely, this means that Y11=o\'^^^\'^ ~ ^ 
forO< fi<n, that A'^^vi^^''^^ = A^^^) if = i^, and that A'^'^vi^^'-^^ = ^ z/ (standard 
linear algebraic theory guarantees the existence of such a basis). Then if the non-zero vector 
B satisfies W^Va"* 7^ 0, we also have that 5° 7^ 0. 

Remark 3.0.2. The "±" sign above distinguishes whether or not most of the eigenvalue of A^'^ 
are positive or negative. For example, if all of the eigenvalues of A^'^ are positive, then Cu = 
Ayu{y) + ■ ■ ■ , while if they are all negative, then Cu = —Ayu{y) + ■ ■ ■ (and similarly for the other 
two cases). 

Proof. Let's consider the first case, in which all of the eigenvalues have the same (non-zero) sign. 
Then by standard linear algebra, since A^^" is symmetric and positive definite (perhaps after mul- 
tiplying it by —1), there exists an invertible "change-of-basis" matrix M^^ such that 

(3.0.24) MJ^A'^^Mf^" = 

where I'^^ '= diag(l, 1, ■ ■ ■ , 1) is the (n + 1) x (n + 1) identity matrix. In fact, we can choose 

(3.0.25) MJ^ = , v^}f^ (no summation in u), 

^ ' ^J\h^\ " ^ ^' 

where A*^'^^ is the "eigenvalue" of A corresponding to the unit-length covector v'^a^ (i.e., ^^=0 I'^cT^P = 
1) appearing in the statement of the theorem. 

We now make the linear change of variables y^ = MJ^x". Then by the chain rule, = f^^^ = 
M/^. Therefore, 

(3.0.26) A"^ — — ,M = A"^M/M/ — — M = I^'' — —u = Ayu. 

This completes the proof in the first case. 

In the second case, in which n of the eigenvalues of A have the same (non-zero) sign, and the 
remaining one has the opposite (non-zero) sign, the proof is similar. The key difference is that 



because of the eigenvalue of opposite sign, (3.0.24) is replaced with 



(3.0.27) M/A°'^M/ = {m-^Y", 



where (m ^)^'^ =^ diag(— 1, 1, 1, ■ ■ ■ , 1) is the standard (1+n) x (1+n) Minkowski matrix. Therefore, 

(3.0.28) A''^ — —u = A''^MJ'MJ' — —u = {m~^Y'' — —u = UyU. 

This completes the proof in the second case. 
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In the third case, in which n of the eigenvalues of A have the same (non-zero) sign, and the 
remaining one is 0, the proof is similar. The key difference is that because of the zero eigenvalue. 



(3.0.24) is replaced with 



(3.0.29) 



where = diag(0, 1, 1, ■ ■ ■ , 1). 
Therefore, 



(3.0.30) 



A 



al3 



d d 



Furthermore, we have that 
(3.0.31) 



u 



. d d 



u = D 



d d 



-u 



i\2 



B'^—u = MJ'B^—u. 



Thus, using using (3.0.25), we have that 



(3.0.32) 



Example 3.0.1. In the first example from above. 



(3.0.33) A'^'^ = 

To calculate the eigenvalues of A, we first set 



1 -2 
-2 2 



(3.0.34) det(A - XI) = det 

The solutions are 



1 - A -2 
-2 2 - A 



A^ - 3A - 2 = 0. 



□ 



(3.0.35) 



A 



3± V17 



Since the eigenvalues are of opposite sign, the corresponding PDE is hyperbolic. 
Example 3.0.2. In the second example from above. 



(3.0.36) Af"" = 

To calculate the eigenvalues of A, we first set 



-2 -1 
-1 -1 



(3.0.37) det(A - XI) = det 

The solutions are 



-2 -A -1 
-1 -1-A 



A^ + 3A + 1 = 0. 
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(3.0.38) 



A 



-3±y5 



Both of these eigenvalues are negative, and thus the corresponding PDE is eUiptic. 
Example 3.0.3. In the final example from above, 



(3.0.39) Ai"^ - 

To calculate the eigenvalues of ^4, we first set 



1 -1 
-1 1 



(3.0.40) det(^ - A/) = det 

The solutions are 



1 - A -1 
-1 1 - A 



A^ + 2A = 0. 



(3.0.41) A = 0, -2, 

and so the corresponding PDE is parabolic. 
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